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Abstract. We derive the lateral Casimir-Polder force on a ground state atom on 
top of a corrugated surface, up to first order in the corrugation amplitude. Our 
calculation is based on the scattering approach, which takes into account nonspecular 
reflections and polarization mixing for electromagnetic quantum fluctuations impinging 
on real materials. We compare our flrst order exact result with two commonly used 
approximation methods. We show that the proximity force approximation (large 
corrugation wavelengths) overestimates the lateral force, while the pairwise summation 
approach underestimates it due to the non-additivity of dispersion forces. We argue 
that a frequency shift measurement for the dipolar lateral oscillations of cold atoms 
could provide a striking demonstration of nontrivial geometrical effects on the quantum 
vacuum. 
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1. Introduction 

The quantum vacuum electromagnetic field plays an essential role in the interaction 
between ground state atoms yy. These van der Waals/Casimir-Polder interactions also 
take place for an atom near a bulk surface or between two bulk surfaces |2] (Casimir 
effect) [3]. Several methods have been employed to measure the atom-surface Casimir- 
Polder interaction [H El El [71 [H [H [TOl [HI [12] . The Casimir force between bulk surfaces 
has been measured for different geometries over the last years [131 HH [13 [IS [IT] , paving 
the way for promising applications in nanotechnology |18j . 

On the theoretical front, surface geometries are usually analyzed with the help of 
the proximity force approximation (PFA, also known as Derjaguin approximation) |19j . 
In this approach, the interaction energy between curved surfaces is obtained from the 
plane-plane case by a simple average over local distances (see, for instance, [20]). PFA- 
based theories have been employed to compare with experimental results for the normal 
Casimir force between a plane and a sphere [I3l|lll|l7], between crossed cylinders [T5] . 
between a plane and a cylinder [21], as well as for the lateral Casimir force between 
corrugated plates [22] . 

However, PFA only holds when the separation distance is much smaller than the 
length scales characterizing the surface geometries (radius of curvature for spherical 
and cylindrical plates, and corrugation wavelength for corrugated surfaces). Thus, it 
is very important to check the accuracy of PFA predictions by comparing them with 
more rigorous formalisms. Beyond-PFA theories for normal Casimir force for cylindrical 
[23| [2l] and spherical [25] plates, and for the lateral Casimir force [261 [22] have been 
recently reported. Most of these theories, though, cannot be directly compared to 
existing experimental data, either because of the assumption of perfect reflectivity 
[23| [211 [26] (a very poor approximation for the typical separation distances of a few 
hundred nanometers employed in the experiments), or because assumptions about 
the geometry itself (small corrugation amplitudes) [27] are not satisfied by the actual 
experimental conditions [22] . 

In this paper, we propose that the lateral Casimir-Polder force on a ground-state 
atom on top of a corrugated surface provides a striking evidence of a non-trivial {i.e. 
beyond PFA) geometry effect. In fact, we show that a Bose-Einstein condensate (EEC) 
trapped close to a corrugated surface may be used as a local field sensor of the geometry- 
induced reshaping of quantum vacuum fluctuations. The control of the atom-surface 
interaction may be relevant for several applications involving cold atoms near surfaces 
|28j . Deviations from the PFA regime are considerably larger than in the case of the 
lateral force for bulk surfaces, where the geometry effect is averaged over the surfaces. 

We employ the scattering approach [29] and write the interaction energy in terms 
of a non-specular reflection operator that takes into account diffraction and polarization 
mixing. A similar formalism was developed previously to analyze the roughness 
correction to the normal Casimir force [30] and to the lateral Casimir force between 
corrugated bulk plates [27J. The corrugation is taken into account as a perturbation 
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of the plane geometry, and the interaction potential is computed up to first order in 
the corrugation amplitude a. Hence our results only hold when a is the smallest 
length scale in the problem: a -C za, ^c, but arbitrary relative values for the atom- 
surface distance za (measured with respect to some arbitrary reference plane) and the 
corrugation wavelength Ac are allowed. Our results are valid for any real material, not 
necessarily a perfect conductor, and exact to first order in perturbation theory in powers 
of the corrugation amplitude. 

We recover the PEA results when taking za ^ Ac- In this limit, the effect is 
dominated by specular reflections, since the the surface is nearly plane in the scale 
of the separation distance. The resulting PFA potential is trivially connected to the 
plane case by taking the local atom-surface distance. We also compare our results with 
those obtained from the pairwise summation (PWS) approach developed in [32]. Since 
the Casimir- Polder force is actually not additive [3], Ref. [32] corrects the result of 
the pairwise integration by a pre-factor "calibrated" by comparing the pairwise and 
exact results for the perfectly reflecting, planar geometry in the retarded Casimir- 
Polder regime {za ^ A^, where A^ is the typical wavelength characterizing the optical 
response of the atom). As could have been expected, the resulting lateral force coincides 
with the more rigorous results presented in this paper when considering nearly plane 
corrugated surfaces, za <^ Ac. In this limit, both our exact result and the approximate 
PWS approach agree with the approximate PFA approach. Outside the PFA regime, 
however, PFA overestimates and PWS underestimates the lateral effect, the discrepancy 
increasing as the corrugation wavelength gets shorter. 

In the present paper, we focus on the comparison between our results and the 
PFA and PWS approaches. For simplicity, we take a perfectly-reflecting surface for the 
numerical examples, leaving the more general discussion to a future publication [31]. A 
short summary of our results, including the discussion of real materials, was reported 
by us elsewhere [35] . 

2. Casimir-Polder potential for a corrugated surface 

The corrugated surface is described by the profile function h{r) (with r = {x,y)) giving 
the local height with respect to the reference plane z = 0. The potential energy of 
a ground-state atom located at some position = {xA,yA, za) near the surface is 
written as 

U{Ra) = U'^'\za) + U'^'\Ra). (1) 

The zeroth-order potential U^^^za), which contains only specular reflection amplitudes, 
is the potential energy for an atom above a plane plate at z = 0. We make a perturbative 
expansion in powers of h and use the scattering approach to find the first-order correction 
U^^\ The first-order scattering by the corrugated surface is described by the matrix 
elements of the first-order reflection operator: (k,p|7?.'-^^|k',p') = i?py (k, k')iJ(k — k'), 
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where H{k.) is the Fourier transform of the profile function h{r). We find 

To first order in the corrugation amphtude this expression is exact, and depends on 
the response function (yf(k, za) which can be exphcitly written in terms of the first-order 
nonspecular reflection coefficients [30] -Rpp'(k, k') as 

g{k,ZA) = ^ / ^ / 7^-^ak',k'-k(^A,0, (3) 

c'eo Jo zTi J (27r)^ 

p,p 

where a{i^) is the dynamic atomic polarizability of the atom (whose ground state is 
assumed to be spherically symmetric), and eo is the vacuum permittivity. Here k and 
k' are the lateral wave vectors of the input and output fields, and p and p' are their 
polarizations (TE for transverse electric and TM for transverse magnetic). The unit 
polarization vectors for outgoing fields are e^E(k) = z x k and e^j^(k) = e^^ ^ 
with K = k + k^z, and = sgn{u) ^ lj"^ / c"^ — k"^. Similar expressions hold for the 
polarization vectors (k) for incoming fields, with the replacement k^ — >■ —kz- The 
roundtrip propagation of the field between the surface and the atom is contained in the 
factor exp[— (k + k')za], with k = y^^^/c^ + fc^. 

Using that the atomic polarizability along the imaginary frequency axis is real and 
positive, one can prove that the response function g{]<., za) is explicitly real. Moreover, 
as expected from symmetry, the response function does not depend on the direction of 
the corrugation profile: (?(k, z^) = g{k,ZA)- 



3. Proximity force approximation 

The response for very smooth surfaces is given by the limit — > in Eq. ([3]). 
According to (HI), (7(0, za) is given entirely in terms of the specular limit of the first-order 
nonspecular coefficients, which quite generally satisfy [30] 

i?pp/(k, k) = 2Krp(k)5py, (5) 

with rp(k) = rp{k) representing the standard Fresnel coefficients for a plane surface. 
This result has a simple interpretation: in the specular limit, the first-order correction 
is just the additional roundtrip phase 2ikzh associated to the surface displacement h. 
After replacing ([5]) into (jl]), we are able to connect the resulting expression for g{0, za) 
with the (zeroth-order) potential for a plane surface: 
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It follows from the "proximity force theorem" that the potential for a very smooth 
surface can be obtained from the planar case by merely taking the local atom-surface 
distance: 

This "proximity force approximation" holds when g{k, za) may be replaced by (7(0, za) 
for all Fourier components of the surface profile function H{k) contributing appreciably 
in Eq. ([2]). Deviations of the PFA from the exact result (for a given Fourier component) 
are measured by the ratio 

_ g{k,ZA) , s 

Ppfa - -77^ r- (oj 

4. Pair wise summation approximation 

The pairwise summation (PWS) technique is an approximation different from PFA that 
computes the Casimir force between bulk surfaces or the Casimir-Polder force between 
an atom and a bulk surface by addition of the atom-atom interaction potential, which is 
U{ri2) = —(23/47r)a^j (0)0:^2(0) ri2 for two ground-state atoms separated by a distance 
ri2 much larger than the static atomic polarizabilities (retarded regime). The PWS 
method has been discussed in [3T] for the normal and lateral Casimir forces between 
corrugated, perfectly reflecting bulk plates, and has been used to calculate lateral CP 
[32] and vdW forces |33] between an atom and a rough or corrugated bulk plate. 

The pairwise lateral Casimir-Polder energy for an atom above a perfectly reflecting 
surface is given by 

U,MRa) = -Cjdr ^^ [(,_,^).^\^_^^p]7 ' (9) 

where C = 15/ica(0)/87r^eo is a "calibration" factor (sometimes improperly called 
"normalization" factor) chosen such that the PWS result coincides with the exact known 
result for a planar surface. To first order in the corrugation amplitude, the lateral PWS 
potential can be written in a form similar to our exact result ([2]) 

ui^si^A) = J ^,e'''-'-gPws{k,ZA)H{k), (10) 
where the PWS response function is 

,pws(k,.^) =-Cj rfrp^^ = -^;^^PWs(fc.A), 

[Z] = e~^(l + Z + ZV3). (11) 



PWS 



Again, due to symmetry, the PWS response function depends only on the modulus of 
k. Deviations of the PWS from the exact result (for a given Fourier component) can be 
measured by the ratio 
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Clearly, in the limit of long corrugation wavelengths kzA <^ 1, both PFA and PWS 
coincide with the exact result. However, as the corrugation wavelength decreases, both 
approximations depart from the exact result given in Eq.(l2]). 

The use of a "calibration" factor C is a major uncontrolled approximation of the 
PWS technique, as it is chosen such that the PWS energy coincides with the exact 
one for a particular distance (in the above, for the retarded Casimir-Polder regime) 
and for a particular model for the surfaces (above, perfect conductors). Clearly, 
such a uniquely fixed calibration factor cannot account for the different geometry and 
materials dependency of dispersion forces in various distances regimes. In view of 
the non-additivity of Casimir forces, the pairwise summation technique is therefore an 
uncontrolled approximation only valid for very dilute media for restricted length scales. 



5. Sinusoidal corrugation 

In this Section we will concentrate on sinusoidal corrugated surfaces, and later on we 
will consider a grooved surface where non-trivial (beyond-PFA) effects of geometry are 
even more impressive. 

For a uni-axial sinusoidal corrugation h{x) = Hq cos{kcx) we find from ([2]) 

U^^\xa,za) = hog{kc,ZA) cos{kcXA)- (13) 

Since g is negative (see below), this potential brings the atom to one of the surface 
crests. For simplicity, we only consider the case of perfect reflectors (see Refs. [35] and 
[3l] for a more general discussion of real materials). In this case, the first-order reflection 
coefficients have a very simple form [30] 

i?TE,TE(k, k') = -2k'C ; i?TE,TM(k, k') = -2^S/c; 

Cf,-' 2kk' 9/^2 
i?TM,TE(k, k') = -2^S ; i?TM,TM(k, k') = + -^C, (14) 

with C = cos{(f) — (f)') and S = sin{(f) — cf)') representing respectively the cosine and the 
sine of the angle between k and k'. These expressions allow us to derive some simple 
analytical results for the exact response function g{k,ZA)- In the non-retarded van der 
Waals regime {za ^ A^) 

hG{kcZA) 



g{kc,ZA) = - ——^ ^ / di a«) 



G{Z) = Z^[2K2{Z) + ZK;{Z)l (15) 

where, as usual, the vdW energy depends on the whole frequency spectrum of dynamic 
atomic polarizability along the imaginary frequency axis. K2 and are the modified 
Bessel functions of second and third order, respectively [36] • In the retarded Casimir- 
Polder regime {za ^ A^) 

3/ica;(0) 

F{Z) =e-^(l + Z+ 16^745 + ZV45). (16) 



g{kc,ZA) = - — ^F{kcZA) 
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As usual, the CP expression depends only on the static atomic polarizability. Both 
analytical results are consistent with the proximity force theorem ([6]) and deviations 
from PEA are given by ppfa = G{kcZA)/ 12 for the vdW regime, and ppfa = FikcZA) 
for the CP regime. Eq. (fT6!) has a greater practical interest because the assumption of 
perfect reflectivity is better suited at longer separation distances. In this Casimir-Polder 
regime, deviations from PWS are given by ppws = -^(^c'2A)/-Ppws(^c-2yi)- 

In Eig. 1 we plot the deviations of the proximity force approximations and of 
the pairwise summation approximation from the exact result. Eor PEA we plot the 
analytical curves for the vdW and CP regime, while for PWS we only plot the curve for 
the CP regime. Since the calibration constant C for the PWS has been calibrated for the 
Casimir-Polder regime, the pairwise summation prediction for the van der Waals regime 
makes no sense, and for this reason we do not plot ppws in the vdW regime. Erom the 
figure we conclude that, for the sinusoidal corrugation, PEA always overestimates the 
lateral force, while PWS always underestimates it. To give a number on the deviation 
from the PEA and PWS, let us consider a separation za = 2/im, well within the 
CP regime, and a corrugation wavelength Ac = 3.5/im {k^ZA ~ 3.55). In this case, 
Ppfa ~ 30% and ppws ~ 115%, which means that the PEA largely overestimates the 
magnitude of the lateral effect, while the PWS underestimates it by as much as 15% for 
this not too large value of kcZA- 

6. Grooved corrugation 

Larger deviations from PEA can be obtained for a surface with periodical grooves (Eig. 
2). In this case the potential (fT3|) has to be replaced by a sum over Eourier components 
a„ of the corrugation profile (assumed to be even for simplicity) h{x) = J2n '^n cos{nkcx), 



oo 



U'^'\xa,za) = J2 



anCosinkcXA) g{nkc,ZA), 



(17) 



n=0 



1.5 - 




vdW regime 
CP regime 
, CP regime 



0.5 - 




k z. 

c A 



Figure 1. Comparison between the first order, exact result with the proximity force 
approximation and the pairwise summation approach for the lateral Casimir-Polder 
interaction between an atom and a perfectly conducting, sinusoidal corrugated surface. 
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■4 ► 

K 

Figure 2. Periodically grooved plate for measuring large corrections to the proximity 
force and pairwise summation approximations. A single cold atom or an elongated 
Bose-Einstein condensate (as in the figure) is trapped with an external asymmetric 
harmonic potential above the plateau region. The dipole oscillation frequency along 
the X direction is shifted due to the lateral Casimir-Polder force. 



where the Fourier coefficients 

{all — for n = , , 

/_^Nn+l 2aK l-cos(nWA.) fo^ 72 > 1 ^ ' 

In Fig. 3 we compare the exact, PFA and PWS lateral Casimir-Polder potentials for 
various values of kcZA- The PFA potential has the same profile as the corrugation itself 
for all values of kcZA-, while the exact and PWS potentials coincide with the PFA version 
only for kcZA <^ 1- The PFA potential is recovered when the response function g{nkc, za) 
may be replaced by the limiting expression g{0,ZA) for all values of n significantly 
contributing to the profile. Otherwise, as za increases, multiplication by g{nkc, za) 
renders the contribution of higher orders comparatively smaller. When kcZA ^ 1, the 
exponential decrease of g implies that the ffist order n = 1 term dominates the sum (ITTll 
(apart from the irrelevant n = term that is x independent and does not contribute to 
the lateral force), and both the exact and PWS potentials are approximately sinusoidal. 
In this limit, the amplitude of the effective sinusoidal exact potential is 

ho = ^ (1 - cos(M/2)) g{k,, Za) ; Ma > 1] (19) 

and for the PWS potential a smiliar equation holds, with g replaced by (^pws- For given 
values of za-, a, and Ac, the amplitude is maximized for s ~ 0.74Ac, which corresponds 
to a rather narrow plateau. It seems more realistic to take, say, s = Ac/2, as we do in 
all our numerical calculations in this paper. 

When the atom is located above one plateau, the PFA predicts that the lateral 
Casimir-Polder force should vanish, since the energy is thus unchanged in a small lateral 
displacement. A non vanishing force appearing when the atom is moved above the 
plateau thus clearly signals a deviation from the PFA. On the other hand, deviations 
from the PWS are not as impressive as for the PFA, but should nevertheless be 
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x(fim) X ([Am) 



Figure 3. Comparison between the exact, proximity approximation, and pairwise 
approximation for the lateral Casimir-Polder potential between an atom and a perfectly 
conducting, periodically grooved surface for different values of kcZA- The groove period 
is Ac = 4/im, the groove width s = Ac/2, and the depth a = 250nm. 



observable for not too large values of kcZA- For example, for kcZA = 3.55, PWS 
underestimates the amplitude of the exact potential by 15%, similarly as for the 
sinusoidal corrugation case. 

7. BEC as a sensor of lateral Casimir-Polder forces 

Dipolar oscillations of the center-of-mass of a Bose- Einstein condensate trapped close to 
a surface have been proposed as a highly sensitive probe of quantum vacuum fluctuations 
[37] . This idea has been successfully implemented to measure the normal component of 
the Casimir-Polder force between rubidium atoms above a dielectric (fused silica) surface 
in the 5 — 10/im separation range |Tl], and to measure non-equilibrium temperature 
corrections to the Casimir-Polder force [I2]. The reported relative frequency shift 
sensitivity in these experiments was 7 = 10~^ — 10^^. 

For the lateral Casimir-Polder force we have in mind a setup where the long axis of 
an elongated BEC would be parallel to the corrugations, while the lateral (along the x 
direction) dipolar oscillation of the center-of-mass would be monitored as a function of 
time (see Fig. 2). When a BEC trapped in an external asymmetric harmonic potential 
of frequencies tu^, = cu^ Uy is set to oscillate perpendicular to the grooves, the lateral 
Casimir-Polder force renormalizes the center-of-mass (CM) dipolar oscillation frequency 
in the x direction, 

2 2 1/", / N d^U^^\x,z) 
^CM,. = ujl + -Jdr no(r) (20) 

where m is the atomic mass, and no(r) = ^^siR"^ — {x'^ + 2;^))^/^ is the two-dimensional 
Thomas-Fermi density for a BEC of Thomas-Fermi radius R. The relative frequency 
shift 7 = {ucu,x - uJx)/uJx is given by 

1 ^k'^ °° 
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X / dp p (1 — p^Y''^ / dO g{nkc, zqm + pR sin 9) cos{nkcpRcos9), 
Jo Jo 

where zcu is the distance from the CM of the BEC to the plate. The relative frequency 
shift 7o for a single atom oscillating above the corrugated surface is obtained in the limit 
R <^ -2cm, Ac, namely 



7o 



(22) 



n=l 



In Fig. 4 we plot the relative frequency shift 7 due to the lateral Casimir-Polder 
force, and compare the exact result with PFA and PWS. Note that the case i? = 
corresponds to the single atom frequency shift 70. As noted before, PFA predicts no 
frequency shift since the potential is locally flat above one plateau. Indeed, 7 is very 
small for kcZciA < 1- However, as fcc-^CM increases, 7 develops a peak and then decays 
exponentially for kcZcM ^ 1- PWS underestimates the lateral force, resulting in a 
frequency shift smaller than the exact one. For the separation chosen in the figure 
(^CM = 2/im), the frequency shift is well above the reported experimental sensitivity 
10~^ — 10^^ [TTII12]. However, the results reported in Fig. 4 correspond to the ideal case 
of a perfectly reflecting plate. As could be expected, the frequency shift is much lower 
for real bulk plates. However, the resulting shifts in those cases are still strong enough 
to allow for a measurement of the lateral Casimir-Polder force with a BEC trapped at 
distances below zqm = Spm [351 Elj- Also, as shown in the figure, the effect of finite 
size of the BEC is to increase the frequency shift as the Thomas-Fermi radius R grows. 



8. Conclusion 

The lateral Casimir-Polder force offers the possibility of testing non-trivial geometrical 
effects of QED. We have calculated the first order, exact lateral force on a ground 




Figure 4. Relative frequency sliift due to the lateral Casimir-Polder force on a 
rubidium BEC (m = 1.45 x IQ-^^kg, a(0)/eo = 47.3 x IQ-^m^) oscillating at a 
distance zcm = 2/xm above a perfectly conducting, periodically grooved surface. The 
exact result, the proximity approximation, and the paiwise summation approximation 
are depicted as a function of kcZcu for different values of the Thomas-Fermi BEC radius 
R. The unperturbed frequency is w/27r — 229Hz, the groove period is Ac = 4/im, the 
groove width s = Ac/2, and the depth a — 250nm. 
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state atom above a corrugated surface and compared it to the widely used proximity 
force and pairwise summation approximations. Our scattering approach aUows one to 
consider real materials and arbitrary atom-surface separations, as long as the corrugation 
amplitude remains the smallest length scale in the problem. On the other hand, the 
PWS approach needs to be "calibrated" for a specific separation distance range since 
the interaction is not additive. Even when comparing in a specific distance range (the 
Casimir-Polder long distance hmit), we have found large deviations from PWS, and even 
larger deviations from PFA. Such deviations could be measured using a Bose-Einstein 
condensate as a surface field sensor. 
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